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Abstract 

We investigate a new "renormalization invariant analytic formulation" of calcu- 
lations in quantum chromodynamics, where the renormalization group summation is 
correlated with the analyticity with respect to the square of the transferred momen- 
tum Q^. The expressions for the invariant charge and matrix elements are then mod- 
ified such that the unphysical singularities of the ghost pole type do not appear at 
all, being by construction compensated by additional nonperturbative contributions. 
Using the new scheme, we show that the results of calculations for a number of physical 
processes are stable with respect to higher-loop effects and the choice of the renormal- 
ization prescription. 

Having in mind applications of the new formulation to inelastic lepton-nucleon 
scattering processes, we analyze the corresponding structure functions starting from 
the general principles of the theory expressed by the Jost-Lehmann-Dyson integral 
representation. We use a nonstandard scaling variable that leads to modified moments 
of the structure functions possessing Kallen-Lehmann analytic properties with respect 
to Q^. We find the relation between these "modified analytic moments" and the 
operator product expansion. 



Take care of the Principles, and the 
Principles shall take care of you. 



Scientific achievements of Nikolai Nikolaevich Bogoliubov are ctiaracterized by a unique 
combination of determination in solving concrete scientific problems and a liigh level of 
matliematical culture. He could find the shortest path to a physical result using most 
general principles of the theory. 

The renormalization-invariant analytic approach to quantum chromodynamics exposed 
here and its most recent applications are based on the works /0, |, ^ by Bogoliubov 
with his closest collaborators. A characteristic feature of these investigations is their strong 
relation with the fundamental quantum physics principles. 
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1 Introduction 



An intrinsic ingredient of modern quantum field theory (QFT) is the renormahzation group 
(RG) method proposed in the mid-fifties [0, |[. The role of this method is particularly 
important in the cases where the interaction is not weak, for example, in quantum chromo- 
dynamics (QCD). Hardly any hadronic process investigated in the QCD framework can be 
analyzed without using the renormahzation group. It is well known that directly solving 
the RG equation for the invariant charge leads to unphysical singularities, for example, to 
the ghost pole in the one-loop approximation. Taking next loop corrections into account 
does not alter the essence, and leads only to additional branch cuts. The existence of such 
singularities contradicts the general principles of local QFT. 

As early as in the late-fifties, N. N. Bogoliubov and collaborators proposed a resolution 
of this problem in the context of quantum electrodynamics (QED) by unifying the RG 
method with the requirement of analyticity with respect to Q^, which in turn followed from 
the known Kallen-Lehmann representation expressing the basic principles of local QFT 
[see Eq. (13) below]. 

The invariant QED charge Q;(Q^) (also referred to as the "invariant, or running coupling 
constant''^) is proportional to the transverse amplitude of the full photon propagator, which 
satisfies the spectral Kallen-Lehmann representation corresponding to the analyticity in the 
complex plane cut along the negative pari^ of the real axis. According to 0, the analytic 
invariant charge can be reconstructed via the Kallen-Lehmann representation, in which the 
relevant spectral density is defined as the imaginary part of the invariant charge determined 
by the RG method in the Euclidean region and analytically continued to the domain where 
ReQ2 < 0. The explicit one-loop (and implicit two-loop) expression obtained in for the 
analytic coupling in QED has the following important properties: 

- the ghost pole is absent; 

- as a function of a, this expression has an essential singularity in the neighborhood of 
a = of the form exp{—3n/a); 

- for real positive a, it admits an expansion in powers of a that coincides with the 
perturbative expansion; 

- it has a finite ultraviolet {UV) limit equal to 3tt, which is independent of the experi- 
mental value a ~ 1/137. 

In P, 0, the idea to combine the renormahzation invariance and the Q^-analyticity in 
QCD led to uncovering new important properties of the analytic coupling. These properties 
include the existence of an infrared fixed point of aa.niQ'^), which proves to be universal 
in the sense that its value a = ^tt / I3q is already determined by the one- loop contribution 
(i.e., remains unchanged by the multiloop corrections and is therefore scheme-invariant). 
It is also independent of the experimentally determined QCD parameter A, and the set 
of curves aan(Q^/A^) corresponding to different values of A is a bundle with the common 
point aan(O) = 47r//3o. Thus, the analytic approach leads to essential modifications of the 
infrared (IR) behavior of the perturbative invariant coupling. We give the approximate for- 
mulas that are useful in the two-loop approximation and also discuss some phenomenological 
applications of the analytic approach.0 

^In view of semantical absurdity of the last term, we use the expression invariant coupling function or 
invariant coupling. 

^We use the notation = —q^, hence the Euclidean region corresponds to positive Q^. 
•^The works [8-16] are devoted to the development and applications of the analytic approach. 
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This work can be conventionally divided into three parts. In the first one (Sec. 2), 
which is a review of our publications over the last two years, the analytic invariant approach 
is formulated in general and is explained in detail in application to the analytic coupling 
"constant." In the second part, which is also a review (Sec. 3), we formulate the "analytic 
perturbation" theory for physical quantities expressed through the two-point objects of the 
type of the Adler D{Q'^) function, whose properties can be related to the Kallen-Lehmann 
representation; we also discuss there the problems of scheme and loop dependence. 

In the third part (Sec. 4), we finally consider the structure functions (formfactors) 
parametrizing the inelastic lepton-hadron scattering cross- sect ion. To relate them to ana- 
lytic functions of Q^, we start with the Jost-Lehmann-Dyson integral representation. Using 
the results of Bogoliubov, Vladimirov, and Tavkhelidze 0], we adduce the arguments in 
favor of the introduction of a special scaling variable such that the moments of the structure 
functions with respect to this variable admit a Kallen-Lehmann representation. This allows 
us to apply the analyticization procedure to these moments. We also consider the relation 
of the analytic moments with the operator product expansion. 



2 An invariant analytic formulation of QCD 

In this section, we formulate the method of constructing the analytic invariant charge and 
consider its main properties. 



2.1 The renormalization group and analyticity 

We start with two remarks. It is known that the invariant QCD charge O-siQ"^) is defined 
via the product of propagators and the special vertex functions, which gives rise to the 
problem of whether the spectral representation can be used for this product. This problem 
was studied in [0, where it was shown that the invariant coupling can be written in the 
form of a spectral integral. In the general case, in addition, the evolution of as{Q'^) is related 
to the "running" gauge parameter. For simplicity, we use the standard MS-scheme, where 
the gauge does not affect the invariant charge.0 

We write the spectral representation for the invariant coupling a((5^) = as((5^)/(47r) as 

aan(Q ) = - / da -. (2.1) 

In the perturbation theory summed up in accordance with the renormalization group, the 
spectral density p(cr, a) decreases as 1/ log^ a, which allows us to write the spectral represen- 
tation without subtractions. 

In the leading logarithmic approximation, the invariant coupling has the form 

aa)fn2) = = \ {2 2) 

^ 1 + a/3o log(QV/i2) /3o log(gVA2) ' ^ • ^ 

where /5o = 11 — 2//3 is the one- loop /5-function coefficient with / active quarks and the 
QCD scaling parameter is A = //exp[— l/2a^/5o]. The corresponding spectral density reads 

"^A similar situation occurs in the MOM-scheme in the transverse gauge or in the MOM-scheme when 
applying a special renormalization p8[. 
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Figure 1: The behavior of the one-loop analytic coupling aan(Q^)- a for A = 200 MeV, b for 
A = 400 MeV. The curves c and d correspond to perturbation theory for the same value 
of A. 



as 

a?l3oTT 1 TT 



Inserting this into spectral integral gives the one-loop analytic coupling function 

A2 



Po 



+ 



log(QVA2) A2-Q2 



(2.4) 



The first term on the right-hand side preserves the standard UV-behavior of the invariant 
coupling. The second term, which comes from the spectral representation and enforces the 
proper analytic properties, compensates the ghost pole at = A^ and is essentially nonper- 
turbative (see the general discussion of this point in ||19|). This term gives no contribution to 
the Taylor series expansion. Thus, the causality and spectrality principles expressed in the 
form of Q^-analyticity, send us the message that perturbation theory is not the whole story. 
The requirement of proper analytic properties leads to the appearance of contributions given 
by powers of that cannot be seen in the original perturbative expansion. We note also 
that unlike in electrodynamics, the asymptotic freedom property in QCD has the effect that 
such nonperturbative contributions show up in the effective coupling function already in the 
domain of low energies and momentum transfers reachable in realistic experiments, rather 
than at unrealistically high energies. 

Thus, synthesis of the renormalization-group invariance and analyticity leads to the an- 
alytic invariant charge without the logarithmic pole and with a finite IR value Q aan(O) = 
47r//3o — 1.396. This limiting value is independent of the experimental information related 
to the normalization point a = a(/i^) or to the parameter A; it is instead determined only by 
the /3-function coefficient related to the general group structure of the Lagrangian. Figure 1 
shows a bundle of curves aaniQ"^) corresponding to different values of A and also the standard 
solutions corresponding to the same A. 

^For numerical estimates at small Q^, we use the number of active quarks / = 3. 
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0,10 




Figure 2: The graph of the one-loop (3-function. 



The graph of the one-loop /5-function illustrating the existence of an infrared fixed point 
in the analytic approach is shown in Fig. 2. The horizontal axis is the parameter j3oa and the 
vertical axis is the function —j3{a). We note that in the one- loop case, one has the symmetry 
with respect to the point PoO' = 1/2, which is broken when taking higher orders into account. 

We now proceed to the two-loop case. The corresponding /3-function reads as 



P{a) = — /3oa^(l + bia), 



hi 



Pi 



(3i = 102 



38/ 



/3o' 3 
Integrating the renormalization group equation, we obtain the transcendental relation 

1 , , / 1 \ 



(2.5) 



Po log X 



a{x) 



hi log 1 + 



hiaix) 



(2.6) 



that can be solved in terms of the Lambert function ||20|, |2T 

0.04 



0.03 



0.02 



0.01 - 



0.00 




Figure 3: The spectral densities: the exact one-loop (a); the exact two-loop (b); the ap- 
proximate two-loop expression (the iterative solution) (c); the exact three-loop function (d). 
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The spectral density obtained from this expression is shown in Fig. 3 (curve b) . It proves 
to be very close to the spectral density corresponding to the explicit iterative solution of 
Eq. O. 

which is useful in the subsequent analysis. 

Solution (|2.7| ) corresponds to the spectral function 



1 + 



rlT\ J.H gi + J^ „ 111 

I [L) = TT + jjiarccos . , Bi = —r . 

'{Bi + L)' + 7r2 /^o 

Its graph is given in Fig. 3 (curve c), where we also show the one- loop (curve a) and the 
three-loop (curve d) results. The three-loop p*-^-* shown in Fig. 3 is obtained in the MS-scheme 
from the exact integral of the RG-equation with the three-loop coefficient 

2857 _ 5033 325 3863 ^ 

^ 2 18 54 6 

As can be seen from Fig. 3, the behavior of spectral densities is stabilized starting with 
the two-loop level; as shown in what follows, moreover, the areas below each of these curves 
are the same, which corresponds to the universality of aan(O). 




Figure 4: Stability of the analytic invariant charge with respect to higher-loop corrections. 
We use the normalization at the T-lepton mass Oan 

(M^) = 0.34 for f = 3; 1 is the one-loop 
approximation, 2 is the two- and three-loop approximations. 
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To obtain a^rHQ'^)^ have to insert spectral density (|2.8| ) in Eq. (|2.1|) . The resulting 
integral cannot be evaluated explicitly.^ The proper analytic properties are reconstructed by 
not only eliminating the pole, but also by subtracting the unphysical branch cut < < 
A^exp(— caused by the double-logarithm dependence in ( p.7| ). 

The numerical calculation results for / = 3 and for the normalization at the point 
aan(M^) = 0.34 are shown in Fig. 4, where we also give the one-loop curve (the corresponding 
values of A are given in Table 1). The three-loop MS-curve is practically identical with the 
two-loop one, with the accuracy of the order 1 %. Thus, in contrast with perturbation 
theory, analyticity leads to an essential stabilization of the invariant charge behavior in the 
IR region. Recalling the asymptotic freedom property, we obtain stability in all the Euclidean 
domain < < oo. 

We note here that the universal behavior of the analytic coupling function is not a 
consequence of the particular two-loop formula ( ^.7| ). The same conclusion remains valid 
when using the exact solution (|2.6|) . Thus, the IR stability of the analytic charge is an 
internal property of the method and is ensured by contributions that are not analytic in as- 
This approach does not introduce any additional parameters into the theory; it operates 
only with the scaling parameter A or with a certain normalization point. 



2.2 Subtraction of unphysical singularities 

The analytic expression for the invariant coupling was obtained using spectral represen- 
tation (|2.1| ) that guarantees the proper analytic properties in the complex plane and 
effectively amounts to subtracting the unphysical singularities (the pole and the cuts). It it 
useful to explicitly separate these terms. 

We consider the complex plane of 2; = Q'^ / h?. The method of subtracting the singularities 
allows us to obtain an explicit expression for the analytic coupling in the one-loop case. 
Indeed, the expression (3Qa'^^\z) = l/log2; has an unphysical pole at z = 1 with the residue 
res [Poa^^\z),z = 1] = 1, whose elimination amounts to adding the term 1/(1 — z), such 
that the expression satisfying the proper analytic properties has the form given in 



z-plane 




Figure 5: Two-loop singularities in the complex plane (z = Q"^ / h?). 



In what follows, we explicitly give the corresponding approximate formulas. 
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In the two-loop case, we first consider ( p.7|) , which in addition to having the ghost pole at 
z = 1 with the residue res [Pod^'^\z), z = 1] = 1/2, has an unphysical cut along the positive 
part of the real axis < 2; < exp(— (see Fig. 5). The subtraction is effected by the pole 
term 

/3oA4oie(^) = (2.10) 



and by the integral 



(3oAd^,{z) = - (2.11) 

TT JO a — z 

TlBi 



[log(a) + B, log(-l - \og{a)/B,)Y + n^Bf 

that eliminates the unphysical branch cut. As the result, the analytic invariant charge can 
be written as 

/3oainH^) = f3od^'\z) + /?oA41(z) + f3oAd^,{z). (2.12) 

This form is convenient because the analytic coupling is represented as a sum of the standard 
expression and of the additional terms of a nonperturbative nature. Their contribution can 
be represented as an expansion in powers of A^/Q^ [see Eq. ( p.22| ) below]. 

2.3 Universality of aan(O) 

The universal value at = is formed by the contribution of the pole term Adp^Y^{z = 
0) = l/(2/5o) and the contribution of (ICTI ) that can be represented as 



dx 1 



PoJo (x + l-logx)2 + 7r2 2po 



The total contribution leads to the universal expression aan(O) = 1/Po- 

In the above approach to approximating the original two-loop coupling, the residue at 
the pole (which is the leading unphysical singularity) is independent of the two-loop f3- 
function coefficient, and it may thus seem that precisely this fact makes aan(O) independent 
of higher-loop corrections. As we have noted, however, there is a different reason behind the 
universality of aan(O), which does not reduce to the choice of a particular approximation of 
the original invariant coupling. We now explain this in more detail. The standard asymptotic 
two-loop expression can be obtained by expanding the function 

Poa^'\z) = ] ^ p , ^ , , ,^,1 , (2.13) 

log z + Bi log[l + log(2;/c;)J 



where C is a constant. Expression ( p. 13 ) correctly reproduces the standard UV limit 

1 A iogiog(gVA2) 



Air 

a. = -5- 

Po 



log(QVA2) PI log^(gVA2 



(2.14) 



that is independent of the constant C. At the same time, the residue at the pole now depends 
on the two-loop /3-function coefficient through Bi, 

res[Poa^'\z),z = I] = ^-1-^ , (2.15) 
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and therefore, the same dependence is involved in the corresponding compensating term 

/3oA.«W = ^^^ (2.16) 



whose contribution to aan(O) is equal to 

^"-.(0) = Arrive' 



The contribution to aan(O) of the term compensating the unphysical branch cut is now 
given by the integral 

. (2) , ^ 1 C r°° dx 



/3o Bi Jo [{x + l)C/Bi - log xf + 7r2 
1 B, 



{Bi + cy 



(2.18) 



which together with the pole contribution ( p.l7| ) gives the universal value aan(O) = l//?o that 
is independent of either C or Bi. 

When taking the higher-loop contribution into account for proving the universality of the 
IR limit in the analytic approach, it is convenient to use the complex quantity ( = 1/a. We 
now give simpler arguments based on the expansion of the perturbative charge into a double 
series in powers of log'"(£)/^'^. For aan(O), we can write 

I foo ^ oo k 

aan(O) = - / dLg{L) = + E E Aafc,^(0) , (2.19) 

^ k=l m=0 

where the higher-loop contribution is given by 

1 /•°° log'^fL - ivr) 
A.,,„(0) = -I,„/_dL^fi-^. (2,20) 

Since the integrand in (|2.2CI|) has no singularities in the lower half-plane, we immediately 
obtain Aafc^m(O) = 0, which proves the universality of the infrared fixed point value of the 
analytic charge. 

Thus, the analyticity requirement for the running charge leads to essential modifications 
of perturbation theory in the IR region. The most relevant factor here is the universality 
of the IR limiting value of the analytic coupling function (the invariance with respect to 
higher- loop corrections), which results in that the family of the invariant charge curves 
corresponding to different loop approximations looks as a bundle with the common point at 

= 0. In addition, these curves obviously come closer to each other in the UV region in view 
of the asymptotic freedom property. In our approach, unlike in the standard perturbation 
theory, there emerges a remarkably stable picture of the invariant charge behavior with 
respect to higher corrections. This stability is important for phenomenological applications, 
where the relevant energy interval is of the order of or less than several GeV. 

2.4 Approximate formulas 



The explicit one-loop formula ( p.4|) is very simple, and its use does not lead to any com- 
plications. In the two-loop case, the analytic coupling is written in the form of an integral 
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representation, and it is interesting to find explicit approximate expressions that are conve- 
nient in applications. 

We consider two such formulas. The first expression follows directly from the picture of 
subtracting the unphysical singularities as explained in Sec. 3. Thus, the analytic coupling 
can be represented as 

«an(g') = apT(g') + A«sing(g') , (2.21) 

where apT(Q^) is a perturbative contribution and the term ^OisingiQ"^) has the effect of 
subtracting the unphysical singularities. For the perturbative term taken as in (p.7|) , the term 
eliminating the unphysical singularities can be represented as two terms whose respective 
effects are to subtract the unphysical pole and the branch cut. The term compensating 
the pole has a simple form. For the term compensating the cut, we use the fact that the 
expansion coefficients C^, 

Aa.M^) ^ -yEXo^jCk, (2.22) 

c = r dt ""M-BiKt + i)] 



(t + l-logt)2 + 7r2 



are numerically small and decrease rapidly (Ci = 0.0354, C2 = 0.0079, C3 = 0.0023, . 
Keeping only the first term in the expansion, we obtain a simple interpolation formula 

^anr)rox(Q ) 



/5o llog(g2/A2) + B, log[l + log(g2/A2)/5^] 

A2 A2 1 

-T^^iK (2-23) 



2 g2 _ A2 Q2 

which provides good approximation [] to the two-loop analytic coupling for moderately 
large Q^. In the interval 1 < Q < 1.5 GeV, the accuracy of the approximation is not 
worse than 0.4%, and for large values of the difference between the formulas becomes 
negligible. Thus, expression ( |2.23D is quite acceptable in the domain of moderately large 



Q > 1 GeV. 

In a number of cases, however, it is necessary to deal with smaller values of Q, down 
to Q — 0. Formula (|2.23| ) is no longer applicable to such problems because the term com- 
pensating the branch cut is poorly approximated by power-series expansion ( p.22| ). The 
approximate formula 



Att 

"approxV^ / n 



1 

+ 



.(g2) i-exp[4(g2) 



"approxW J 

.XQ') = In ^ + 5i In iL ^ + 4vr2 



(2.24) 



for the two- loop analytic charge can be used also for Q — 0. Equation ( p.24| ) reproduces the 
UV two-loop asymptotic behavior ( p.l4| ) and the universal limiting value at = 0. This 



expression approximates the exact one for Q > 1 GeV with the accuracy within 1 % and can 
be used for all Q^. 

^The approximate formula for the two-loop correction [a^]an to the physical quantities of the Z?-function 
type can also be found in this way. 
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Table 1: The perturbative and analytic one- and two-loop values of the scaling parameter 
(MeV) for / = 3 versus the normalization point as{M^). 



«s(M2) 


0.30 


0.32 


0.34 


0.36 


0.38 




173 


201 


228 


256 


283 


Ad) 

an 


197 


235 


275 


319 


366 




333 


377 


419 


460 


500 


A(2) 
an 


434 


516 


607 


706 


814 


A(2) 
appr 


423 


500 


582 


671 


777 



For sufficiently large Q^, the analytic coupling function is dominated by its perturbative 
component. Already for Q = M^, however, the nonperturbative contribution becomes es- 
sential. In Table 1, we compare the A parameter values corresponding to the perturbative 
and analytic approaches. The result obtained according to Eq. (|2.23|) reproduces the ex- 
act two-loop calculation with high accuracy and is not given given here. For the two-loop 
perturbative formula, we used expression (|2.7|) , which is most appropriate for our analysis.^ 
The bottom row corresponds to approximate expression ( p. 24 ). 



3 Analytic perturbation theory 

In this section, we briefly review applications of the analytic approach to the analysis of 
several processes. For the physical quantities considered here, we use the analyticization 
procedure of the entire perturbative expression involving higher powers of the invariant 
charge 0. This strategy leads to the so-called analytic perturbation theory (APT). 

We consider the integral characteristics of the invariant charge in the IR region by extract- 
ing the relevant information from the physics of jets, and also from the e^e~-annihilation 
processes into hadrons and the inclusive r-lepton decay. We use this set of data to study 
the dependence of theoretical results on the choice of the renormalization scheme. We show 
that applying the APT allows us to considerably reduce the scheme dependence. This in 
turn means that the three-loop level attained for many processes is practically independent 
of the choice of the scheme. 



3.1 The integral characteristics of as in the IR region 

A distinctive feature of the analytic charge is that it is flnite in the IR region. This property, 
which is sometimes referred to as the coupling "freezing," is often used for phenomenological 
purposes (see, for example, the discussion in |2^). Experimental evidence for the regular IR 
behavior of the QCD charge was ingeniously extracted from physics of jets using the integral 
characteristics 

A{Q) = ^ dkas{k^). (3.1) 
Q Jo 

^We note that using formula (O) as a perturbative one, leads to somewhat greater values of A than 



when working it out from Eq. (2.14) 
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It has been empirically found ^ that A{2 GeV) = 0.52 ± 0.10. 



Table 2: The infrared integral characteristics of as{k^) evaluated in the one- and two-loop 
approximations for normalization at the r-lepton mass. 



«an(M2) 


0.34 


0.36 


0.38 


^l-loop(2 ) 


0.50 


0.52 


0.55 


^2-loop(2 ) 


0.48 


0.50 


0.52 



We normalize at the r-lepton mass. Calculations of A{2 GeV) are given in Table 2. 
It can be seen that the APT approach allows us to uniformly and consistently describe 
the almost-perturbative region of the order of the r-lepton mass and the nonperturbative 
characteristics ( |3.1| ) without introducing any additional parameters. 



3.2 The e+e -annihilation process into hadrons 

We now apply the analytic approach to the analysis of the e"'"e~-annihilation process into 
hadrons. To compare the results with the experimental data, we use the method of so called 
"smearing" of resonances proposed in The analysis of the e"'"e~-annihilation into hadrons 
carried out in ||2^ relied on a certain "optimum" renormalization scheme constructed on the 
base of the principle of minimal sensitivity (PMS) p6[ with the third-order perturbation 
theory used for optimization. Our analysis is not based on any optimization of the scheme 
arbitrariness. Moreover, we show that the scheme dependence in the APT is considerably less 
than in the standard approach, and its predictions have practically no scheme arbitrariness 
in the entire energy range. 

The analyticization procedure can be also applied to observable quantities for which the 
appropriate analytic properties are known. The APT can be applied to an object that has 
numerous applications, namely the Adler D-function 



-Q 



f 



(3.2) 



where n(s) is the correlation function and d{Q'^) is the QCD correction that is expanded in 
the RG perturbation theory as 



d{Q^) = a(Q2)[l + d,a{Q^) + d^a^Q^) + ■ ■ 



(3.3) 



where f\a = agju. 

The D-function is related to the function -R(s) defined as the ratio of the hadron and 
lepton cross-sections for the e"'"e~-annihilation by 



/•c 

Dm = Q' / 

Jo 



ds 



Q 



2^2 



R{s) 



(3.4) 



^In this formula, we allowed ourselves to change the normalization of the coupling constant so as to 
simplify comparing with the previous works on the subject, where, as a rule, the quantity a ~ ols/t^ is used 
as the invariant charge. 
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This also implies the properties of D{Q'^) as an analytic function in the Q^-plane cut along 
the negative semi-axis. We define the spectral density p'^^(o") through the discontinuity 
of ( |3.3|) on this cut, 

p^ff(a) = p(i)(a) + t/ip(')((T) + d2P^^\a) + ■■■. (3.5) 



The expression p''^-'(cr) is the spectral function of the invariant charge and p^^\a) in 
corresponds to the kill power of the effective coupling. Thus, the analytic expression for the 
QCD correction to the D-function is written as 

d^MQ^) = 4pt(Q') + did^lUi.Q'') + d25tU{Q'') + ■ • • , (3.6) 

where the first term 5apt(Q^) coincides with the analytic invariant charge. The subsequent 
terms do not reduce to powers of the analytic coupling; thus, the APT method leads to 



non-power-series expansions. Properties of such expansions were analyzed in |T2 . 

We define the QCD correction r(s) to the function R{s) in the same manner as for the 
D-function in (13.21), and we use the relations 



dm = Q' r r^K^), r{s) = ^d{-z) , (3.7) 

where the integration contour in the last expression is in the analyticity domain of the 
integrand and bypasses the cut along the real semi-axis. 

We take the quark thresholds into account by using the approximate formula proposed 
in H, 

R{s) = 3j2Q}e{s- 4m}) T{vf) [1 + g{vf) rf{s)] , (3.8) 
/ 

where Vf and the functions T{v) and g{v) are given by 



4:71 \ n 3 + f/7r 3 
9iv) = - 



(3.9) 



2v 4 \2 4tx 

In the APT, the correction r/(s) is expressed through the effective spectral density as 

r,(.) = 1 r-pf{a\ (3.10) 

where p'f^ip) is defined in terms of the discontinuity of dflQ"^) on the physical branch cut. 
The corresponding three-loop contribution is written as 

df{Q') = a/(Q2)[i + d«a/(g2) + d^pajiQ% (3.11) 

where the MS-scheme coefficients are equal to 

= 1.986 - 0.115/, 
df = 18.244 - 4.216 / + 0.086 /2 + rffs^'^* , 



_ 1.2395 (E'fQf 



14 



It is hardly possible to use perturbative expressions for a direct description of the exper- 
imentally observed quantity R{s), because of the threshold singularities of the form (a<j/f 



We use the "smearing" method proposed in |24], which does nevertheless allow us to compare 



the results with the experiment. The idea of this approach consists in replacing the quantity 
R{s) defined through the correlation function 11 as 

R{s) = ^^[U{s + te)-U{s-ze)], (3.12) 

with the quantity 

= ^ [n(s + zA) - n(s - iA)] (3.13) 

for some finite A. For the values of s near the threshold, quantity ( |3.12| ) is very sensitive 
to the threshold singularities, in the vicinity of which the perturbative expansion looses 
its applicability. Stepping away from the real axis into the complex plane by a finite 
distance A, as in (|3.13| ), we can expect that it would be possible to describe ( |3.13|) using an 
appropriate perturbative approximation. 

The "experimental" curve corresponding to ( p.l3| ) can be found if we use the dispersion 



relation for the correlator 11 (g^) to write Eq. ( p.l3| ) as 



^a(3) = - r ds'^^^—-. (3.14) 
71 Jo (s — s')^ + A^ 



The corresponding "experimental" curves were found in [22| for some values of A, whose 
estimates were made in [^. We use these curves for comparing with our results. 

We note that a direct use of perturbation theory for describing Ra{s) is again impossi- 
ble. Indeed, the i?-ratio in (|3.14|) parametrized using the invariant charge with unphysical 



singularities leads to a divergence of the integral in ( p.l4| ). Thus, even though the use of the 
"smeared" quantity ( 3.14|) allows us to bypass the complication with the threshold singular- 



ities, there arises a problem related to the behavior of the running charge in the IR region. 
We can avoid this complication using the APT. 

For A = 3 GeV^, Fig. 6 shows the corresponding experimental curve and the curve found 
in ||22| from the PMS-optimization of the third-order perturbative expansion. The same 
figure gives also the result of our calculation through the third order.0 For the scaling 
parameter in the analytic approach, we took the value Aan = 870 MeV (/ = 3) obtained 
from the analysis of the semileptonic r-decay in the APT framework. For the quark masses, 
we took the values that are close to the constituent ones (cf. [0), wi„ = rrid = 250 MeV, 
m, = 400 MeV, = 1.35 GeV, = 4.75 GeV and rrit = 17 A GeV. 

New "experimental" data for the D-function were obtained recently We give in 



Fig. 7 the corresponding curves and also the result of our calculation. Figure 7 shows that 
good fit of the experimental data can be achieved already in the first order of APT. The same 
conclusion, as we have noted, is valid for R/^ in the entire energy range IQ. We note here 
that loop stability is not observed in the standard approach using the PMS-optimization. 
Moreover, the whole "trick" is based here on higher approximations. Thus, the situation 
regarding the absence or the presence of the infrared fixed point that can emerge in scheme 
optimization of the perturbative expansion depends in an essential way on the quantity under 
consideration (i.e., is defined by the coefficients of the perturbative expansion) p9 . 



""^^As shown in the calculation of -Ra in the analytic approach leads to good fit of the experimental 
curve already in the first order. 
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Figure 6: The quantity R/^{q^) corresponding to the parameter value A = 3GeV^. The 
figure shows the experimental curve, the result of the PMS-optimization of the third-order 



perturbative expansion obtained in J2^/ , and the result of the analytic approach through the 
third order. 
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Figure 7: D- function. 
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3.3 The dependence on the renormahzation scheme 

Inevitable termination of the PT series, i.e., the approximation of a physical quantity by 
one of its partial sums, leads to the known problem of the dependence of the results on the 
renormahzation prescription. Thus, the partial sum of the PT series used in approximating 
a physical quantity bears a dependence on the choice of the renormahzation scheme, which is 
the source of theoretical ambiguity in describing experimental data. In QCD, such ambiguity 
is the greater the smaller are the energy parameters characteristic of the process. To solve 
the stability problem of the results obtained, it is by far not enough to investigate only 
loop stability within a certain renormahzation scheme; one should also consider the scheme 
stability of the results. 

We discuss the scheme arbitrariness arising in the APT in the example of the i?-ratio 
for the e"*"e~-annihilation process into hadrons. We consider a class of MS-like schemes and 
compare our results with those obtained in the perturbative analysis (see, for example, []30| ). 

In passing from one renormahzation scheme to another, the coupling constant transforms 

as 

a = a(l + Via + t>2a^ + ■■■). (3.15) 

We limit ourselves here to the three- loop level of the D-function achieved at present, with 
the QCD corrections taken in the approximation where 

d = a{l + dia + d2a^) , (3.16) 

with the running charge determined as a solution of the renormahzation group equation with 
the three-loop /3-function 

da 

13(a) = /i^TT^ = -6a^(l + ha + 62^^) , (3.17) 

where 



62 



^ ^ 33-2/ 153-19/ 

6 66 -4/ ' ^ ' 

MS _ 77139 - 15099/ + 325/2 



288(33 - 2/) 



The three-loop /3-function coefficient 62 and the expansion coefficients di and d2 in (|3.16| ) 
depend on the choice of the renormahzation scheme. Under scheme transformation ( |3.15| ), 
they change as 

b'2 = 62 - '^1 - hvi + V2 , 

d[ = di-vi, (3.19) 
d2 = d2 - 2{di - vi)vi - V2 . 



Thus, every term in representation ( |3.16| ) undergoes a transformation, and we thus obtain 
the new function 

d' = a'{l + d[a' + d'^a'^), (3.20) 

where the coupling a' is evaluated with the new /3-function, with the three-loop coefficient 62 
replaced by the primed one 62- 
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Recalling the transformation law of the scaling parameter 

A' = Aexp{vi/b) 
and Eqs. ( |3.19| ), we find two scheme invariants 

b Q2 

pi = -log— - c/i, p2 = h2 + d2 - hidi - d\. (3.21) 

We normalize the momentum scale at A^. In arbitrary scheme, the invariant charge is 
then determined from the equation 

=rfp-rfi + $(a,62), (3.22) 

where 




^ , , , 1 , , 1 + 6ia /■'^ dx . . 

$(a, 6,) = - - log + h2 1 jTTb^xjiTTh^Th^y ^ ^ ^ 

Although there are no general arguments to prefer a certain renormalization scheme from 
the start, we nevertheless can define a class of "natural" schemes, which look reasonable at 



the three-loop level that we consider. The relevant criterion was proposed in ||32[. One 
should restrict oneself to the schemes where the cancellations between different terms in the 
second scheme invariant ( |3.21| ) are not too large. Quantitatively, this criterion can be related 
to the cancellation index 

C = T^i\b2\ + \d2\+dl + \di\bi). (3.24) 

One should of course keep in mind the conventions involved in these considerations, in 
particular as regards the minimal value of the cancellation index. 

Given a certain maximum value of the cancellation index Cmax? we can investigate sta- 
bility of the results obtained by taking different schemes with the index C < Cmax- As Cmax, 
we take the index corresponding to the optimal PMS-scheme. We then have a relatively 
small class of "admissible" schemes bounded by the maximal index CpMs- 

For R{s), the cancellation index Cr is evaluated using the known coefficients ri and r2 of 
the perturbative expansion of the correction r = a(l -|- ria + r2a^). For the PMS-scheme, it 
is CpMs — 2. To demonstrate the scheme arbitrariness arising here, we choose two schemes 
from this class. The first one is the H scheme with the parameters rf^^ = —3.2 and 62^^ = 
(the 't Hooft scheme), and the second is the MS-scheme corresponding to the parameters 
^(MS) _ -j^ ^2^^'' = 4.47. These schemes are close to each other and to the boundary 

cancellation index Ch — Cy^ ~ C'pms — 2. 

Figure 8 shows the QCD correction r(s) as a function of y/s/Ajj^ evaluated in perturba- 
tion theory and in the analytic approach for two renormalization schemes H and MS with 
approximately the same cancellation indices Cr ~ 2. As can be seen from the figure, the 
analytic approach allows us to drastically reduce the scheme arbitrariness. 

Essential reduction of the scheme dependence in the APT also takes place for other pro- 
cesses, for example the inclusive r-decay [jll|, and in the Bjorken and Gross-Llewellyn Smith 



sum rules for the inelastic lepton-hadron scattering |T3| , |T^ . In the analytic approach, there- 
fore, the three-loop level reached presently for a number of physical processes is practically 
invariant with respect to the choice of the renormalization prescription. 
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Figure 8: The graph of r{s) calculated in perturbation theory (PT) and in the analytic 
approach (AA) for two renormalization schemes H and MS with the same cancellation 
index Cji^2. 

3.4 Inclusive r-lepton decay 

The inclusive r-decay (see Fig. 9 for the corresponding diagram) allows one to perform a 
low-energy test of QCD. The r-lepton mass Mr = 1777l;o;26 MeV [Q, on the one hand, is 
sufficiently large to allow the hadronic decay modes, but on the other hand, is small in the 
chromodynamics scale, where it is in the low-energy domain. Theoretical description of the 
inclusive r-decay is in principle possible without any model assumptions, which is important 
for reliably determining the low-energy value of as(M^) from experimental data. The main 
quantity to be studied is the i?T--ratio 

Rr = + (3.25) 

r[r Z^e(7)J 

which in the modern experiments can be measures with the accuracy of several per cent. 
The starting point of the theoretical analysis is the expression 




Rr = 2i \l + T7^\R{s\ (3.26) 



where R{s) is defined by the imaginary part of the hadron correlator 

n(s)= E IK/[n,gy(s) + n„,,A(s)]. (3.27) 

q=d,s 

Here Vuq are the Kobayashi-Maskawa matrix elements. In the massless case considered here, 
the vector and axial-vector hadron correlators, H-uqy and Iluq^A respectively, coincide, and 
the function R{s) is equal to the ratio R{s) for the e"'"e~-annihilation process into hadrons. 



19 




Figure 9: The inclusive T-lepton decay diagram. 

The standard analysis of the r-decay immediately faces a difficulty in applying the original 
formula ( p.26 ), because the parametrization of the function R{s) by perturbative as with the 



unphysical singularities leads to singularities in the integrand. The way out proposed in 
consists in the following. Integral ( p.26|) is represented as a combination of integrals along the 
sides of the cuts in the s complex plane (see Fig. 10). By the Cauchy theorem, this integral 
is then "transformed" into the integral along the contour |s| = M^. After the integration by 
parts, we are left with the contour representation of Rr involving the D-function, 

R^ = ^<f —il-z)\l + z)D{M^z). (3.28) 
2m J\z\=i z 



The transition from the original expression ( 3. 26 ) to contour representation ( 3.28 ) is 



based on certain analytic properties of the correlator, which are violated in the standard 
analysis. Thus, the proper analytic properties ensuring the analytic approach are important 
for the consistency of the inclusive r-decay description. 

We describe this process in the APT 0. We single out the strong-interaction contribu- 
tion A,- to the -Rr-ratio 

i?, = (1 + A,), (3.29) 

where is a known factor including electroweak corrections. 
We express A,- through the effective spectral function as 

di da di r^'r da ( V/ \ 



A. = ^ / -P^"(-) - -1-77^1 + 77^ P^"(-)- (3.30) 



TT Jo a ' ' n Jo a \ J \ M^j' 

Because of the universality property, the integral in the first term can be expressed through 
aan(O). The spectral function in the two- loop approximation has the form 



where the spectral density of the invariant charge g{a) is defined in ( p.8|) and the functions 



R{L) and /(L) are given by ( |2.9|) . Inserting ( |3.31|) in (|3.30| ) allows us to evaluate the strong- 



interaction contribution A,- in terms of the scale parameter A. 



Using the experimental value Rr = 3.633 ± 0.031 |3|, we obtain a{M^) = 0.400 ± 0.026 



and the corresponding value of the scaling parameter A^^^ = 935 ± 159 MeV. These values 
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Figure 10: Transition to the contour representation for R^-. 



are larger than those obtained in PT using the contour representation The reason hes 
in the fact that the nonperturbative corrections characteristic of the analytic approach give 
a negative contribution to |T^. Thus, to obtain the same value in PT and in 

the analytic approach, the "perturbative component" contribution of the latter should be 
increased by increasing A. The inclusive r-decay was analyzed at the three-loop APT level 
36[|. The corresponding value of A turned out to be smaller, A^'^^ = 871 ± 155 MeV. 



m 



The scheme stability of this analysis was also shown in It should be noted that the 



quantity Aan is very sensitive to the experimental value of Rr- Thus, using R^. = 3.559±0.035 
|37| , we obtain A^^^ = 640 ± 127 MeV, which corresponds to a considerably smaller invariant 
charge at the mass (see Table 1). 



4 The analytic approach in inelastic lepton-hadron 
scattering 

In this section, we give a theoretical foundation of a possible application of our anal5d;ic de- 
scription to inelastic lepton-hadron scattering processes. The key point of our construction — 
the analytic properties of the structure function moments with respect to — requires a 
certain modification of the standard formalism, in particular, the change of the standard 
Bjorken moments M„(Q^) with the modified moments AiniQ"^) with respect to a new scal- 
ing variable that takes kinematic mass dependence into account. We start with the Jost- 
Lehmann integral representation (see, for example, §55 of f^) for the Fourier image of the 
corresponding matrix element. 
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4.1 The Jost— Lehmann representation 



The structure functions of the inelastic lepton-hadron scattering depend on two arguments, 
and the corresponding representations that accumulate the fundamental properties of the 
theory (such as relativistic invariance, spectrality, and causality) have a more complicated 
form in our analysis than in representations for functions of one variable. Two such rep- 
resentations are known in the literature. We use the 4-dimensional integral representation 
proposed by Jost and Lehmann |]38[ for the so-called symmetric case.|^ Applications of this 
representation to automodel asymptotic structure functions were considered by Bogoliubov, 
Vladimirov, and Tavkhelidze |^ , some of whose results and notation we use in what follows. 
The proof of the Jost-Lehmann representation is based on the most general properties of 
the theory, such as covariance, Hermiticity, spectrality, and causality (see 0; some math- 
ematical problems related to the Jost-Lehmann-Dyson representation are also considered 
11). 



m 
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For definiteness, we speak about the inelastic scattering of charged leptons (electrons, 
muons) on nucleons, i.e., we consider the process £ + — *• £ -|- hadrons. In the lowest order 
in the electromagnetic coupling constant (one-photon exchange), this process is shown in 
Fig. 11, which also explains our notation. In the unpolarized case, the cross-section of the 
process is defined by the hadronic tensor 



W^^{q, P) = -^Y^ j dx exp{iq ■ x) (^P, a 



Ml 



1 Jy 



(4.1) 



constructed of the commutators of the currents, with the sum taken over the nucleon polar- 
izations. 




Figure 11: The deep inelastic lepton-hadron scattering diagram in the one-photon exchange 
approximation . 

Relativistic invariance and the electromagnetic current conservation lead to the parametriza- 
tion of tensor (^?l| ) in terms of two structure functions Wi and W2, 



(4.2) 



+ -THlPf^- \Pu - ^qu] W2iq, P) 



^^A more general case was considered by Dyson [|9|, and similar representations are therefore often called 
the Jost-Lehmann-Dyson representations. 
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where M = v -P^ is the nucleon mass. 

We now hst the main properties of the functions w following from the general principles 
of local QFT: 

- covariance property means that the functions w depend on two scalar arguments, which 
we choose as u = P ■ q and = — g^, 

w{q,P) = W{u,Q')- 

- spectrality property is written as 

W{u,Q^) = for ^=x>l, 

where we used the dimensionless Bjorken variable, which in the physical domain of the 
process for (g + P)^ > is kinematically restricted by the interval < a; < 1; 

- the structure function parametrizes the scattering cross-section and is real (the reality 
property), 

W{u,Q^)=W*iu,Q'); 

- Hermiticity of the current operator leads to the (anti-) symmetry property 

Wi-u,Q') = -Wiu,Q'y, 

- the vanishing of the commutator of currents at space-like intervals because of the local 
commutativity of currents gives the causality condition 

/ exp{-iqz)W{q, P) = for < g. 

For the function W^(z/, Q"^) satisfying all these conditions, there exists a real moderately 
growing distribution ip{u, A^) such that the Jost-Lehmann integral representation holds; in 
the nucleon rest frame, this can be written as ||^ 

W{u,Q^) =e{qo) JdudX^6[ql-{Mu-qf-X^]^{u,X^), (4.3) 

with the function ip{u, A^) supported on the set 



p = |u|<i, > M\i- ^1- p^y. 

For the process under consideration, the physical values of u and are positive. We, 
thus, can neglect the factor e{qo) = e^u) and keep the same notation for W{h',Q'^). Taking 
into account that the weight function ip{u,X'^) = il){p,X^) is radial-symmetric, as follows 
from covariance |Q , we write the Jost-Lehmann representation for W in the covariant form, 

Wiu,Q^) = dpp^ dX^ / dz (4.4) 

JO Jx^ . J-1 



mill 



X (5(g2 + mV + A^ - 2zp^iy^ + M^Q^ )^{p, A^) , 

where 

>L. = M\l-Jl^^f. (4.5) 
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4.2 Analytic moments of the structure functions 

As follows from representation (|4.4| ), a natural scaling variable is given by 



^ 1 / Q2(Q2 + 4M2) ^ / g2 + 
^~2V Z/2 + M2Q2 ~ Q2 + 4]\^2a.2 ' ^ 

which accumulates the root structure determined by the (5-function argument. At the same 
time, in the physical region of the process, the s variable changes in the same way as the 
Bjorken variable x, i.e., from zero to one. The variable s bears a dependence on the mass of 
the target (the nucleon) and is different from both the Bjorken variable and the Nachtmann 
variable pf 



1 + ^1 + xHMyq^ 

that is sometimes used in the kinematical account of mass effects of inelastic scattering pro- 
cesses. However, only the s variable leads to the moments that have the analytic properties 
in that we need. 

To establish these properties, we integrate over z in (|4.4|) as 



1 poo 

2 



W{iy,Q') = , / dpp dy (4.7) 



2 

mill 



X 6 



V^(P,A^). 



'1 + 4MVQ2 

We now define the modified s-moments of the structure functions (cf. 



^'^(^^) = n^,MVn^,in-.y. I dss^-'Wiu, Q^) . (4.8) 



;i + 4MVQ2) 
Inserting W{u, Q^) as given by ( [4.7|) , we obtain 



n Jo Jo [Q^ + a)"- 

where a = + M^p^ and a^m = 2M^{1 - p^). 
Introducing the weight function 

mn{a) = - f' dpp"^'e{a-a^i^)tlj{p,a-M^p^), (4.9) 



n Jo 

we obtain the representation for the s-moments 

M„m ^ f (4,10) 

which implies the analyticity of AiniQ"^) in the complex plane cut along the negative 
semi-axis, i.e., the Kallen-Lehmann type analyticity. 
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In [Q, the Deser-Gilbert-Sudarshan integral representation |^ was used to arrive at a 
similar statement regarding the analyticity of the Kallen-Lehmann type for the x-moments. 
However the status of this representation in QFT is less clear, since it cannot be obtained 
starting with only the basic principles of the theory. 

The relation between moments ([4.8|) and the standard Bjorken moments 



Mn{Q^)= f'dxx''-^W{u,Q^) (4.11) 
Jo 



can be expressed by 



M(Q^) = y nk + {n + l)/2] 

^^^^ r[(n + l)/2] k\ 

M(Q^) = f nfc + (n + l)/2] 

'^^^^ nin + l)/2]to k\ 
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In the asymptotic domain corresponding to large values of the transferred momentum Q 
where power corrections of the form 1/(Q^)" can be neglected, the x-, s-, and ^-moments are 
identical. Outside the asymptotic domain, on the other hand, where it comes to studying 
the contribution of higher twists, the difference between these definitions of moments must 
be taken into account. 



4.3 Dispersion relation and the operator product expansion 

To establish the relation with the operator product expansion, we start with the Jost- 
Lehmann representation and obtain a dispersion relation for the forward Compton scattering 
amplitude with respect to the new variable ( [4.6|) . We write the matrix element of the process 
corresponding to representation ( [4.4|) as 



1 /•! roo /•! 

T(z/,g2) = _ dpp^ dz (4.12) 

TT JO . J-1 

mm 

In the complex plane, the function T(z/, Q^) has a branch cut along the positive part 
of the real axis starting at v^^^ defined by the condition 

g2 + + A 



vl,^ + = niin 



2zp 

Recalling ( [4.5|) and the range of the integration variable z in (^4.12|) , we can simplify this 
to ^ 

^rrw^-ny '^^^^^^'V^' . (4.13) 

{p} 2p 

which leads to z/^j^ = {Q^/2f. Thus, the sought dispersion relation has the form 

1 roo dip' 

TM') = -l, 2 I ■ W{v,,Q'). (4.14) 
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We note that in terms of the Bjorken variable x, relation ( 4.14 ) is represented as 



This expression determines simple properties of the amplitude T{x, Q^) in the complex x- 
plane and is convenient in the operator product expansion. 

In considering consequences of the Jost-Lehmann representation, as noted above, the 
natural scaling variable is given by s. In this case, there arises a similar structure of the 
dispersion integral 



n Jo si 1 — {Si/s)^ 



(4.16) 



4s^ 

The identity between the structures of the dispersion relations with respect to the vari- 



ables X and s allows us to establish the relation of analytic moments (|4.8| ) to the operator 
product expansions of currents used in finding the Q^-evolution of the structure functions 
of the moments. The moments in Eq. ( [4.11|) correspond to the case where only the Lorentz 
structures of the form P^^ . . . P^^ are taken into account in matrix elements of the opera- 
tor {P\0^^,,,^J\P). Then the application of the operator product expansion for the Compton 
amplitude leads to the expansion in powers of {q-P)/Q'^, i.e., to the expansion in the inverse 
powers of x. A similar expansion in the inverse powers of x can also be done in dispersion 
integral ( 4.15|) . The coefficients are then determined by the x-moments. Comparing the 



two power series gives the sought relation between the x-moments and the operator product 
expansion. 

In the general case, the symmetric matrix element (-P|0^^...^„ |P) contains the Lorentz 
structures given by {P^^ . . . P^„}, M'^g^^^^{P^^ . . . -P^^.j}, etc. The moments with respect 
to the ^ variable correspond to choosing the operator basis where the expansion goes over 
traceless tensors, i.e., such that the contraction of g^^^. with (-P|0^^...^^ |P) vanishes for any 
two indices. It is then obvious that the Lorentz structure of the matrix element (-P|0^^...^^ |P) 
is fixed unambiguously. 

Dispersion representation ( [4.16 ) allows us to expand the Compton amplitude in the 



inverse powers of s. If the operator basis is chosen such that an arbitrary contraction of the 
tensor {P\0 ^^,„^^\P) with the nucleon momentum P^^ vanishes, then the operator product 
expansion leads to a power series for the forward Compton scattering amplitude with the 
expansion parameter q^qu[P^iPu — 9 fifiP"^) / {q'^Y , which corresponds to expanding dispersion 
integral ( 4.16|) in powers of 1/s^. We thus arrive at the relation between the analytic s- 



moments and the operator product expansion. We stress that the orthogonality requirement 
of the symmetric tensor (P|0^^...^,JP) to the nucleon momentum P^. determines its Lorentz 
structure unambiguously. 



5 Conclusions 

We considered the analytic formulation of QCD, where the analyticized RG-solutions for 
the invariant coupling functions, the Green's functions, and the matrix elements are free 

^^We note that in using other scahng variables, for example, the Nachtmann one, this structure can be 
destroyed. 
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of unphysical singularities. An important property of this formulation that we found is the 
stability of the analytic invariant charge with respect to higher-loop correction in all of the 
range. The key point here is the existence of the universal limiting value Q;an(0) = ^t^/ I^q that 
is invariant with respect to multiloop corrections. This constant is independent of the Aqcd 
parameter and is determined only by the general symmetry properties of the Lagrangian. 
Therefore, the family of curves aan(Q^/A^) for different values of the A parameter is a bundle 
with the common point Q;an(0) = 47r/ /3o (this picture is independent of the number of loops). 

The invariant analytic formulation essentially modifies the behavior of as_^{x) in the IR 
region by making it stable with respect to higher-loop corrections. The two-loop approxi- 
mation differs from the one-loop one by no more than ~ 10 % in the small-Q^ domain, and 
the three-loop approximation differs from the two- loop one by only ~ 1 %. This is radically 
different from the situation encountered in the standard renormalization-group PT, which is 
characterized by strong instability with respect to the next loop corrections in the domain 
of small ~ A^. We note also that maintaining the proper analytic properties with re- 
spect to is essential for a self-consistent definition of the effective coupling function in the 
time-like region |[T^ . In describing the concrete processes, for example the inclusive r-lepton 



decay, a consistent analysis is possible only provided the above analytic properties hold. 

There are at least two possibilities to describe physical quantities in the new approach 
framework. The simplest one consists in replacing as{x) — as_^{x) in the explicit expres- 
sions for the observables "processed" by the RG method, or more precisely, for the related 
quantities defined in the space-like region of the variable. 

We take another possibility, however. For the quantities similar to the Adler D{Q^)- 
function that are represented by the PT power series, according to a special convention, 
the analyticization procedure is applied to each power of asiQ"^) separately. This leads to 
a new non-power-series expansion, in which the powers of asiQ"^) are replaced with new 
nonsingular functions of We call this algorithm, which was first proposed in 0, 

the APT. Applying this algorithm to analyze the amplitudes of the processes like the e~^e~- 
annihilation into hadrons and the inclusive r-decay, and also of the sum rules for the inelastic 
lepton-hadron scattering, we see that in addition to possessing loop stability, the APT 
results are much less sensitive to the choice of the renormalization scheme than in the 
standard approach. In other words, the three-loop APT level practically insures both the 
loop saturation and the scheme invariance of the relevant physical quantities in the entire 
energy or momentum range. 

It appears that by accounting for the additional information about the proper analytic 
properties, the first terms of the APT non-power-series expansion already give sufficiently 
good approximation to the sum of the whole series. We recall here the analogy with summing 
up the perturbative expansions with the additional information on the behavior of the remote 
PT series terms taken into account [^. In that case, it also turned out that the expression 
for the approximated function given by the first several terms of the loop expansion was 
practically unchanged by higher corrections. 

In this work, we considered also the structure functions of the inelastic lepton-hadron 
scattering, which are more complicated objects than the two-point functions, which are in 
one way or another related with the Kallen-Lehmann representation. For these functions, 
the general quantum field theory principles, including covariance, Hermiticity, spectrality, 
and causality, are expressed by the Jost-Lehmann-Dyson integral representation. In using 
the analytic approach to define the Q^-^volution, it was convenient to introduce the moments 
■MniQ'^) of the structure functions corresponding to the special scaling variable (|4.6|). It is 
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these moments, rather than the Bjorken or Nachtmann ones, that exhibits simple analytic 
properties with respect to Q^. In this work, we found the relation of the new analytic 
moments A^„(Q^) to the operator product expansion, where the tensor structure of the 
matrix elements of operators with respect to the nucleon states must be fixed according to 
the condition that they be orthogonal to the nucleon momentum. 
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